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Prediction and control of the dynamics of complex networks is a central problem in network 
science. Structural and dynamical similarities of different real networks suggest that some universal 
laws might accurately describe the dynamics of these networks, albeit the nature and common 
origin of such laws remain elusive. Here we show that the causal network representing the large- 
scale structure of spacetime in our accelerating universe is a power-law graph with strong clustering, 
similar to many complex networks such as the Internet, social, or biological networks. We prove 
that this structural similarity is a consequence of the asymptotic equivalence between the large- 
scale growth dynamics of complex networks and causal networks. This equivalence suggests that 
unexpectedly similar laws govern the dynamics of complex networks and spacetime in the universe, 
with implications to network science and cosmology. 

I. INTRODUCTION 

Physics explains complex phenomena in nature by reducing them to an interplay of simple fundamental laws. This 
very successful tradition seems to experience certain difficulties in application to complex systems in general, and to 
complex networks in particular, where it remains unclear if there exist some unique universal laws explaining a variety 
of structural and dynamical similarities found in many different real networks [TM71. One could potentially remedy 
this situation by identifying a well-understood physical system whose large-scale dynamics would be asymptotically 
identical to the dynamics of complex networks. One could then try to use the extensively studied dynamical laws of 
that physical system to predict and possibly control the dynamics of networks. At the first glance, this programme 
seems to be quite difficult to execute, as there are no indications where to start. Yet we show here that there exists 
a very simple but completely unexpected connection between networks and cosmology. 

In cosmology, de Sitter spacetime plays a central role as the exact solution of Einstein's equations for an empty 
universe, to which our universe asymptotically converges. Here we show that graphs encoding the large-scale causal 
structure of de Sitter spacetime and our universe have structure common to many complex networks fS^lOJ , and that 
the large-scale growth dynamics of these causal graphs and complex networks are asymptotically the same. To show 
this, we describe the causal graphs first. 

The finite speed of light c is a fundamental constant of our physical world, responsible for the non-trivial causal 
structure of the universe [TT]. If in some coordinate system the spatial distance x between two spacetime events 
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FIG. 1: Finite speed of light c, and causal structure of spacetime. In panel (a), a light source located at spatial coordinate 
a; = is switched on at time £ = 0. This event, denoted by L in the figure, is not immediately visible to an observer located at 
distance lo from the light source. The observer does not see any light until time t — xo/c. Since no signal can propagate faster 
than c, the events on the observer's world line, shown by the vertical dashed line, are not causally related to L until the world 
line enters the L's future light cone (yellow color) at t = xq/c. This light cone depicts the set of events that L can causally 
influence. An example is event P located on the observer's world line x = a;o at time t = to > xo/c. The past light cone of P 
(green color) is the set of events that can causally influence P. Events L and P lie within each other's light cones. Panel (b) 
shows a set of points sprinkled into the considered spacetime patch. The red and green links show all causal connections of 
events L and P in the resulting causet. These links form a subset of all the links in the causet (not shown). 



(points in space and time) is larger than ct, where t is the time difference between them, then these two events cannot 
be causally related since no signal can propagate faster than c (Fig. flja)). Causality is fundamental not only in 
physics, but also in fields as disparate as distributed systems [T21[T3] and philosophy [TJ. 

The main physical motivation for quantum gravity is that at the Planck scale {Ip ~ 10"'^^ meters and tp ^ 10^'*^ 
seconds), one expects spacetime not to be continuous but to have a discrete structure |15| . similar to ordinary matter, 
which is not continuous at atomic scales but instead is composed of discrete atoms. The mathematical fact that 
the structure of a relativistic spacetime is almost fully determined by its causal structure alone |16H18| motivates 
the causal set approach to quantum gravity [19J . This approach postulates that spacetime at the Planck scale is a 
discrete causal set, or causet. A causet is a set of elements (Planck-scale "atoms" of spacetime) endowed with causal 
relationships among them. A causet is thus a network in which nodes are spacetime quanta, and links are causal 
relationships between them. To make contact with General Relativity, one expects the theory to give rise to causal sets 
which are constructed by a Poisson process, i.e. by sprinkling points into spacetime uniformly at random, and then 
connecting each pair of points iff they lie within each other's light cones (Fig.IIJb)). According to the theorem in |20| . 
causets constructed by Poisson sprinkling are relativistically invariant, as opposed to regular lattices, for example. 
Therefore we will use Poisson sprinkling here to construct causets corresponding to spacetimes. An important goal in 
causal set quantum gravity (not discussed here) is to identify fundamental physical laws of causet growth consistent 
with Poisson sprinkling onto realistic spacetimes in the classical limit [211 112] ■ 

In 1998 the expansion of our universe was found to be accelerating J23 [ l24 j . Positive vacuum energy, or dark energy, 



corresponding to a positive cosmological constant A in the Einstein equations (Bl), is currently the most plausible 
explanation for this acceleration, even though the origin and nature of dark energy is one of the deepest mysteries 
in contemporary science j25j . Positive A implies that the universe is asymptotically (at late times) described by de 
Sitter spacetime [551 [57]. We first consider the structure of causets sprinkled onto de Sitter spacetime, then quantify 
how different this structure is for the real universe, and finally prove the asymptotic equivalence between the growth 
dynamics of de Sitter causets and complex networks. 

II. RESULTS 
A. Structure of de Sitter causets 

De Sitter spacetime is the solution of Einstein's field equations for an empty universe with positive cosmological 
constant A. The 1 + 1-dimensional de Sitter spacetime (the first '1' stands for the space dimension; the second '1' — 
for time) can be visualized as a one-sheeted 2-dimensional hyperboloid embedded in a flat 3-dimensional Minkowski 
space (Fig. [2la)). The length of horizontal circles in Fig. [2|a), corresponding to the volume of space at a moment of 
time, grows exponentially with time t. Since causet nodes are distributed uniformly over spacetime, their number also 
grows exponentially with time, while as we show below, their degree decays exponentially, resulting in a power-law 
degree distribution in the causet. 

To obtain this result, we consider in Fig. [21 a) a patch of 1 + 1-dimensional de Sitter spacetime between times t — 
(the "big bang") and t = to > (the "current" time), and sprinkle N nodes onto it with uniform density 6. In 
this spacetime the element of length ds (often called the metric because its expression contains the full information 
about the metric tensor) and volume dV (or area, since the spacetime is two-dimensional) are given by the following 



expressions (Section B 1 ) 



ds^ == -dt^ +cosh^tde^, (1) 

dV = cosht dtde, (2) 

where 9 e [0, 27r) is the angular (space) coordinate on the hyperboloid. In view of the last equation and uniform 
sprinkling, implying that the expected number of nodes dN in spacetime volume dV is dN = S dV, the temporal node 
density p{t) at time t e [0, to] is 

, , cosht f f 

Pit) = ^^^ ^ e'-'\ (3) 

smh to 

where the last approximation holds for tg > t ^ 1. 

Since links between two nodes in the causet exist only if the nodes lie within each other's light cones, the expected 
degree k{t) of a node at time coordinate t G [0, to] is proportional to the sum of the volumes of two light cones centered 
at the node: the past light cone cut below at t = 0, and the future light cone cut above at t = to, similar to Fig. [T] 
Denoting these volumes by Vp(t) and Vy(t), and orienting causet links from the future to the past, i.e. from nodes 
with higher t to nodes with lower t, we can write k{t) = ko(t) + ki(t), where ko{t) = SVp{t) and fci(t) = SVf{t) are 




FIG. 2: Mapping between the de Sitter universe and complex networks. Panel (a) shows the 1 + 1-dimensional de Sitter 
spacetime represented by the upper half of the outer one-sheeted hyperboloid in the 3-dimensional Minkowski space XY Z. The 
spacetime coordinates (S,i), shown by the red arrows, cover the whole de Sitter spacetime. The spatial coordinate 6'o of any 
spacetime event, e.g. point P, is its polar angle in the XY plane, while P's temporal coordinate io is the length of the arc lying 
on the hyperboloid and connecting the point to the XY plane where f = 0. At any time i, the spatial slice of the spacetime is a 
circle. This 1-dimensional space expands exponentially with time. Dual to the outer hyperboloid is the inner hyperboloid — the 
hyperbolic 2-dimensional space, i.e. the hyperbolic plane, represented by the upper sheet of a two-sheeted hyperboloid. The 
mapping between the two hyperboloids is shown by the blue arrows. The green shapes show the past light cone of point P in 
the de Sitter spacetime, and the projection of this light cone onto the hyperbolic plane under the mapping. Panel (b) depicts 
the cut of panel (a) by the Y Z plane to further illustrate the mapping, shown also by the blue arrows. The mapping is the 
reflection between the two hyperboloids with respect to the cone shown by the dashed lines. Panel (c) projects the inner 
hyperboloid (the hyperbolic plane) with the P's past light cone (the green shape) onto the XY plane. The red shape is the left 
half of the hyperbolic disc centered at P and having the radius equal to P's time io, which in this representation is P's radial 
coordinate, i.e. the distance between P and the origin of the XY plane. The green and red shapes become indistinguishable 
at large times io as shown in panels (d,e,f) where these shapes are drawn for io = 5, 10, 15 using the exact expressions from 
Section^ Assuming the average degree of fc = 10, these io times correspond to network sizes of approximately 40, 200, and 
2000 nodes. 
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FIG. 3: Degree distribution in the universe. Panel (a) shows the rescaled distribution Q{k,tq) = Sa'^P{k,to) of rescaled 
degrees k = k/{Sa'^) in the universe causet at the present rescaled time tq — to/a — 0.85, where S is the constant node density 
in spacetime, and a — \/3/A. As shown in Section C the rescaled degree distribution does not depend on either S or a, so 
we set them to 5 = 10* and a = 1 for convenience.^ The size A'' of simulated causets can be also set to any value without 
affecting the degree distribution, and this value is A'^ = 10® nodes in the figure. The degree distribution in this simulated 
causet is juxtaposed against the numeric evaluation of the analytical solution for Q{k., tq) shown by the blue dashed line. The 
inset shows this analytic solution for the whole range of node degrees k £ [1, 10^**] in the universe, where S ~ 10^^^ and 
a ~ 5 X lO'^^. Panel (b) shows the same solution for different values of the present rescaled time r, tracing the evolution of the 
degree distribution in the universe in its past and future. All further details are in Sections | A 5|C| 

the expected out- and in-degrees of the node. One way to compute Vp{t) and V/(i) is to calculate the expressions for 
the light cone boundaries in the {t, 9) coordinates, and then integrate the volume form dV within these boundaries. 
An easier way is to switch from cosmological time t to conformal time [26j 7y(i) = arcseccoshi. After this coordinate 
change, the metric becomes conformally flat, i.e. proportional to the metric ds^ = —dt^ + dx^ in the flat Minkowski 
space in Fig. [T] 

ds^ = sec^ rj [-dyf 
dV — sec^ rjdijdO, 



d9^] 



so that the light cone boundaries are straight lines intersecting the coordinate (77, ( 
{t,x) replaced by {r],6). Therefore, the volumes can be easily calculated: 

rn{t) rvW-n' 

Vp{t) = / drj' / (i6' sec^ 77' = In sec 77 (t) = In cosh i « t, 

Jo Jo 

rn{to) /•'?'-'?(<) 



(4) 
(5) 
-axes at 45°, as in Fig. [I] with 
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(7) 



where approximations hold for io ^ ^ ^ li and where we have used 'q{t) = arcseccoshi w 7r/2 — 2e"*. For large times 
t ^ 1, the past volume and consequently the out-degree are negligible compared to the future volume and in-degree, 
which decay exponentially with time i, 



k,{t) = 5Vf{t) « 5e^°- 



(8) 



These results can be generalized to d -I- 1-dimensional de Sitter spacetimes with any d and any curvature K = 
A/3 = 1/a^, where a, the inverse square root of curvature, is also known as the curvature radius of the de Sitter 
hyperboloid, or as its pseudoradius. Generalizing Eqs. (3|8), we can show that the temporal density of nodes and 
their exnecteH in-decrree in this case scale as p"(*~*o) anH pWo— *) with n- = /? = din.. Tn short, we have a, rnmbination 



their expected in-degree in this case scale as e'^ 
of two exponentials, number of nodes 



d/a. In short, we have a combination 



" born at time t and their degrees ^ e " . This combination yields a 
power-law distribution P{k) ^ k^'^ of node degrees k in the causet, where exponent 7 = 1 + a/ /3 = 2. 



B. Structure of the universe and complex networks 



The large-scale causet structure of the universe in the standard model differs from the structure of sparse de Sitter 
causets in many ways, two of which are particularly important. First, the universe is not empty but contains matter. 
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FIG. 4: Degree distribution and clustering in complex networks and de Sitter spacetime. The Internet is the network 
representing economic relations between autonomous systems, extracted from CAIDA's Internet topology measurements [32j . 
The network size is A'^ = 23752 nodes, average degree k = 4.92, and average clustering c = 0.61. Trust is the social network 
of trust relations between people extracted from the Pretty Good Privacy (PGP) data |33]; A'^ = 23797, k = 7.86, c = 0.48. 
Brain is the functional network of the human brain obtained from the fMRI measurements in |34|; N = 23713, k — 6.14, 
c = 0.16. De Sitter is a causal set in the 1 + 1-dimensional de Sitter spacetime; N — 23739, k = 5.65, c — 0.82. Panel (a) 
shows the degree distribution P{k), i.e. the number of nodes N{k) of degree k divided by the total number of nodes N in the 
networks, P{k) — N{k)/N, so that k = ^f, kP{k). Panel (b) shows average clustering of degree-fc nodes c(fc), i.e. the number 
of triangular subgraphs containing nodes of degree k, divided by N{k)k(k — l)/2, so that c = J]]^. c{k)P{k). All further details 
are in Sections I A HA 4| 



Therefore it is only asymptotically de Sitter |26l [27], meaning that only at large times t ^ a, or rescaled times 
T = t/a 3> 1, space in the universe expands asymptotically the same way as in de Sitter spacetime. In a homogeneous 
and isotropic universe, the metric is ds^ = —dt^ + B?{t) d^^, where dfl is the spatial part of the metric, and function 
R{t) is called the scale factor. In de Sitter spacetime, the scale factor is R{t) ^ coshr, while in a flat universe 
containing only matter and dark energy, R{t) ~ sinh ' (3r/2). In both cases, R{t) ^ e^ at large times r ^ 1, but at 
early times t < 1 the scaling is different. In particular, at t — 7> the universe scale factor goes to zero, resulting in 
a real big bang. The second difference is even more important: the product between the square of inverse curvature 
a'^ = 1/ K^ and sprinkling density 5 = l/{lptp) (one causet element per unit Planck 4-volume) is astronomically huge 
in the universe, Sa'^ ~ 10^^'*, compared to Sa''^'^^ ^ 1 in sparse causets with a small average degree. Collectively these 
two differences result in that the present universe causet is also a power-law graph, but with a different exponent 
7 = 3/4 (Fig. ^a)). 

However, the 7 = 2 scaling currently emerges (Fig.pTb)) as a part of a cosmic coincidence known as the "why now?' 
puzzle |28H31j . The matter and dark energy densities happen to be of the same order of magnitude in the universe 
today. This coincidence implies that the current rescaled time tq = to /a is approximately 1. Figure p[b) traces the 
evolution of the degree distribution in the universe in its past and future. In the matter-dominated era with r < 1, the 
degree distribution is a power law with exponent 3/4 up to a soft cut-off that grows with time. Above this soft cut-off, 
the distribution decays sharply. Once we reach times t ^ 1, e.g. today, we enter the dark-energy-dominated era. The 
part of the distribution with exponent 3/4 freezes, while the soft cut-off transforms into a crossover to another power 
law with exponent 2, whose cut-off grows exponentially with time. The crossover point is located at kcr ^ Sa^. Nodes 
of small degrees k < kcr obey the 7 = 3/4 part of the distribution, while high-degree nodes, k > kcr, lie in its 7 = 2 
regime. At the future infinity r — > 00, the distribution becomes a perfect double power law with exponents 3/4 and 
2. 

In short, the main structural property of the causet in the present-day universe is that it is a graph with a power-law 
degree distribution, which currently transitions from the past matter-dominated era (r < 1) with exponent 7 = 3/4 
to the future dark-energy-dominated era (r > 1) with 7 = 2. In many (but not all) complex networks the degree 
distribution is also a power law with 7 close to 2 [5til0j. In Fig. Wa) we show a few paradigmatic examples of 
large-scale technological, social, and biological networks for which reliable data are available, and juxtapose these 
networks against a de Sitter causet. In all the shown networks, the exponent 7 « 2. This does not mean however that 
the networks are the same in all other respects. Degree-dependent clustering, for example (Fig. llFb)), is different in 
different networks, although average clustering is strong in all the networks. Strong clustering is another structural 
property often observed in complex networks: average clustering in random graphs of similar size and average degree 
is lower by orders of magnitude [8HIO] . 



C. Dynamics of de Sitter causets and complex net^vorks 

Is there a connection revealing a mechanism responsible for the emergence of this structural similarity? Remarkably, 
the answer is yes. This mechanism is the optimization of trade-offs between popularity and similarity, shown to 
accurately describe the large-scale structure and dynamics of some complex networks, such as the Internet, social 
trust network, etc. |3S] The following model of growing networks, with all the parameters set to their default values, 
formalizes this optimization in [35j . New nodes n in a modeled network are born one at a time, n = 1, 2, 3, . . ., so 
that n can be called a network time. Each new node is placed uniformly at random on circle S^. That is, the angular 
coordinates 0„ for new nodes n are drawn from the uniform distribution on [0, 27r]. Circle §^ models a similarity 
space. The closer the two nodes on §^, the more similar they are. All other things equal, the older the node, the more 
popular it is, the higher its degree. Therefore birth time n of node n models its popularity. Upon its birth, new node 
no optimizes between popularity and similarity by establishing its fixed number m of connections to m existing nodes 
n < no that have the minimal values of the product nA9, where A6 = tt — Jtt — |6'„ — 6'„„|| is the angular distance 
between nodes n and uq. One dimension of this trade-off optimization strategy is to connect to nodes with smaller 
birth times n (more popular nodes) ; the other dimension is to connect to nodes at smaller angular distances A9 (more 
similar nodes). After placing each node n at radial coordinate r„ — Inn, all nodes are located on a two-dimensional 
plane at polar coordinates {rn,dn)- For each new node uq, the set of nodes minimizing nA9 is identical to the set of 
nodes minimizing x — \ii(nnoA9/2) = r„ -|- r„g -I- ln(A0/2), where x is equal to the hyperbolic distance [36] between 
nodes hq and n if r, tq, and A6 are sufficiently large. One can compute the expected distance from node ng to its 
mth closest node, and find that this distance is equal to ln[7r?TT,no/{2(l — l/rio)}] ~ r„u + ln[7rTO/2] w r„g, where 
approximations hold for large uq. In other words, new node ng is born at a random location on the edge of an 
expanding hyperbolic disc of radius tq = Innoj and connects to asymptotically all the existing nodes lying within 
hyperbolic distance tq from itself. The connectivity perimeter of new node no at time no is thus the hyperbolic disc 
of radius vq centered at node uq. The resulting connection condition x < vq, satisfied by nodes n to which new node 
no connects, can be rewritten as 

AO 
r„ + In — < 0. (9) 

This model yields growing networks with power-law degree distribution P{k) ^ k ' and 7 = 2. The networks in 
the model also have strongest possible clustering, i.e. the largest possible number of triangular subgraphs, for graphs 
with this degree distribution. The model and its extensions describe the large-scale structure and growth dynamics of 
different real networks with a remarkable accuracy |35j . We next show that the described network growth dynamics 
is asymptotically identical to the growth dynamics of de Sitter causets. 

To show this, consider a new spacetime quantum P that has just been born at current time i = to in Fig- [2| That 
is, assume that the whole de Sitter spacetime is sprinkled by nodes with a uniform density, but only nodes between 
t = and t = to are considered to be "alive." We can then model causet growth as moving forward the current time 
boundary t = to one causet element P at a time. By the causet definition, upon its birth, P connects to all nodes 
in its past light cone shown by green. As illustrated in Fig. [2] we then map the upper half of the outer one-sheeted 
hyperboloid representing the half of de Sitter spacetime dS^ with t > 0, to the upper sheet of the dual two-sheeted inner 
hyperboloid, which is the standard hyperboloid representation of the hyperbolic space H^ [37]. This mapping sends a 
point with coordinates {t, 0) in dS^ to the point with coordinate (r, 9) in H^, where r = t. Since in the conformal time 
coordinates the light cone boundaries are straight lines intersecting the (77, 0)-axes at 45°, the coordinates {t, 6) of all 
points in P's past light cone satisfy inequality A9 < Arj = rj{to) — ri{t) = arcsec cosh to — arcsec cosh i « 2 (e~* — e~*°). 
If to 3> t ^ 1, then we can neglect the second term in the last expression, and the coordinates (t„,0„) of existing 
causet nodes n to which new node P connects upon its birth are given by 

A6'<2e-*", (10) 

which is identical to Eq. Q since r„ = t„. In Section [b] we fill in further details of this proof, extend it to any 
dimension and curvature, and show that the considered mapping between de Sitter spacetime and hyperbolic space 
is relativistically invariant. 

In short, past light cones of new nodes, shown by green in Fig. [21 are asymptotically equal (Fig. [2ld-f)) to the 
corresponding hyperbolic discs, shown by red. The green light cone bounds the set of nodes to which node P connects 
as a new causet element. The red hyperbolic disc bounds the set of nodes to which P connects as a new node 
in the hyperbolic network model that accurately describes the growth of real networks. Since these two sets are 
asymptotically the same, we conclude that not only the structure, but also the growth dynamics of complex networks 
and de Sitter causets are asymptotically identical. 



III. DISCUSSION 

Geometrically, this equivalence is due to a simple duality between the two hyperboloids in Fig. [2J The inner 
hyperboloid represents the popularity x similarity hyperbolic geometry of complex networks; the outer hyperboloid is 
the de Sitter spacetime, which is the solution to Einstein's equations for a universe with positive vacuum energy. In 
that sense, Einstein's equations provide an adequate baseline description for the structure and dynamics of complex 
networks, which can be used for predicting network dynamics at the large scale. De Sitter spacetime is homogeneous 
and isotropic, as is the hyperbolic space, but if we take a real network, e.g. the Internet, and map it to this homogeneous 
space, then after the mapping, the node density in the space is non- uniform |38| . In real networks, the space thus 
appears homogeneous only at the largest scale, while at smaller scales there are inhomogeneities and anisotropics, 
similar to the real universe, in which matter introduces spacetime inhomogeneities at smaller scales, and leads to 
non-trivial coupled dynamics of matter density and spacetime curvature, described by the same Einstein equations. 
In view of this analogy, equations similar to Einstein's equations may also apply to complex networks at smaller 
scales. If so, these equations can be used to predict and possibly control the fine-grained dynamics of links and nodes 
in networks. 

Our results may also have important implications for cosmology. In particular, de Sitter causal sets have exactly 
the same graph structure that maximizes network navigability 39J. Translated to asymptotically sparse causal sets, 
does this property imply that the expanding portion of de Sitter spacetime (t > 0) is the spacetime that maximizes 
the probability that two random Planck-scale events have an ancestor in their common past? If it does, then this 
uniqueness of de Sitter spacetime may lead to a different perspective on the cosmic coincidence problem, as well as 
on dark energy, possibly casting the latter as a phenomenon emerging from certain optimization principles encoded 
in the causal network structure. 

The degree distributions in some complex networks deviate from clean power laws, the exponents of these power 
laws vary a lot across different real networks, and so do clustering, correlation, and many other structural properties 
of these networks [SHIP]. Therefore it may seem unlikely that de Sitter causets can model the full spectrum of 
structural diversity observed in complex networks. Focusing on the trust network in Fig. [4] for instance, we have 
already observed that its degree-dependent clustering is quite different from the one in de Sitter causets. Yet, given 
that these causets are asymptotically identical to growing hyperbolic networks, this observation appears as a paradox, 
because the hyperbolic networks were shown to accurately match not only clustering peculiarities, but also a long list 
of other important structural properties of the same trust network, as well as of other networks |35j . The explanation 
of this paradox lies in that the hyperbolic network model has parameters to tune the degree distribution exponent, 
clustering strength, node fitness, and other network properties, while in de Sitter causets, only the number of nodes 
and average degree can be controlled. Do the hyperbolic network parameters have their duals in the de Sitter settings, 
what are the physical meanings of these dual parameters, and do they lead to similar modeling versatility — all these 
questions are open. 

We conclude with the observation that the node density in growing hyperbolic networks with the default parameters 
corresponding to de Sitter causets, is not uniform in the hyperbolic space [3S]. This observation means that these 
networks are not random geometric graphs [40], and that their structure does not exactly reflect the geometry of 
the underlying hyperbolic space. Informally, a random geometric graph is a coarse, discrete representation of a 
smooth geometric space. Our finding that asymptotically the same networks have a uniform node density in de Sitter 
spacetimes dual to hyperbolic spaces, strongly suggests that real networks are random geometric graphs that grow in 
spacetimes similar to the asymptotically de Sitter spacetime of our accelerating universe. 
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Appendix A: Data and Simulations 



Here we describe the network data and methods used in Figs.[3]^4] The considered real networks are paradigmatic 
examples of complex technological, social, and biological networks for which reliable large-scale data are available. 
We note that in all the three considered networks, links represent soft relational data instead of hard-wired network 
diagrams: economic/business relations in the Internet, trust relations between people, causal/correlation relations in 
the brain, and causal relations in causal sets. 

1. Internet data 

The Internet topology in Fig.Hrepresents economic/business relations between autonomous systems or ASs. The AS 
is an organization or an individual owning a part of the Internet infrastructure. The network is extracted from the data 
collected by CAIDA's Archipelago Measurement Infrastructure (Ark) [32], http://www.caida.org/projects/ark/ 
The Ark infrastructure consists of a set of monitors continuously tracing IP-level paths to random destinations 
in the Internet. The union of the paths collected by all the monitors is then aggregated over a certain period 
of time, and each IP address in the collection is mapped to an AS owning this address, using the RouteVicws 



BGP tables http : //www . routeviews . org/( The resulting AS network has a power-law degree distribution with 



exponent 7 = 2.1, and this result is stable over time [4T1|42] and across different measurement methodologies [43 l 144 ) . 



http://www.caida.org/research/topology/topo_comparison/. The data and its further description are available 



at http://www.caida.org/data/active/ipv4_routed_topology_aslinks_dataset.xml. In Fig. |4j the data for 
June 2009 is used. The aggregation period is one month, and the number of monitors is 36. 



Trust data 



The trust network in Fig. |4] represents trust relations between people, extracted from the Pretty-Good- 
Privacy (PGP) data. PGP is a data encryption and decryption computer program that provides cryptographic 



privacy and authentication for data communication [33], http://www.openpgp.org/. In the PGP trust net- 
work, nodes are certificates consisting of public PGP keys and owner information. A directed link in the net- 
work pointing from certificate A to certificate B represents a digital signature by the owner of A, endorsing the 
owner/public key association of B. We use the PGP web of trust data collected and maintained by Jorgen Ced- 
erlof http: //www. lysator. liu. se/~jc/wotsap/wots2/. We post-process the data as follows. The directed network 
is first mapped to its undirected counterpart by taking into account only bi-directional trust links between the certifi- 
cates. The largest connected component is then extracted from this undirected network. This way we strengthen the 
social aspect of the network, since we consider only pairs of users (owners of PGP keys) who have reciprocally signed 
each other's keys. Such filtering increases the probability that the connected users know each other, and makes the 
extracted network a reliable proxy to the underlying social network. The degree distribution in the resulting network 
is a power law with exponent 7 = 2.1, and this result is stable over time |35j , http: //pgp- cs .uu.nl/plot/ . In Fig.BJ 
the data snapshot taken on December 1, 2005, is used. 

3. Brain data 

The brain functional network in Fig. ffl represents causal/correlation relations between small areas in the human 
brain. This network is extracted from the fMRI measurements in |34j . In those experiments, the whole brains of 
different human subjects are split into 36 x 64 x 64 = 147456 adjacent areas called voxels, each voxel of volume 
3 X 3.475 X 3.475 mm^. The subjects are then asked to perform different tasks, during which the magnetic resonance 
activity V{x,t) is recorded at each voxel x at time t. Time t is discrete: 400 recordings are made with the interval of 
2.5 s. Given this data, and denoting by (•) the time average, the correlation coefficient r{x, x') for each pair of voxels 
is then computed: 

^ {V{x,t)V{x',t))-{V{x,t)){V{x',t)) 

y^[{V{x, ty) - {V{X, i))2] [{V{x',t)^) - {Vix', i))2] ■ 

To form a functional network out of this correlation data, two voxels x and x' are considered causally connected if the 
correlation coefficient between them exceeds a certain threshold r^ r{x,x') > Tc- If Tc is too small or too large, then 
the resulting network is fully connected or fully disconnected. There exists however a unique percolation transition 
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FIG. 5: (a) The total versus in-degree distributions in the de Sitter causet of Fig. B] (b) The out-degree distribution in the 
same causet. The solid line shows the Poissonian distribution with the mean X = ko = k/2 = 2.83, where ko is the average 
out-degree in the causet. 



value of Tc corresponding to the onset of the giant component in the network. To find this value, we compute the sizes 
1 5*1 1 r^ and |S'2|r<: of the largest and second largest components 5*1 and ^2 in the network for different values of r^. The 
value of threshold Tc corresponding to the largest values of |9|S'i|/9rc|re and \S2\r^ is then its percolation transition 
value. With threshold re set to this value, the network has a power-law degree distribution with exponent 7 = 2 
and exponential cutoffs, and this result is stable across different human subjects and different types of activity that 
they perform during measurements, see [SUHS] and http://www.caida.org/publications/papers/2012/network_ 
[cosmology/supplemental/. In Fig. I4J a specific dataset is used — Setl4, see the URL — where the subject is at rest. 
The threshold value is Tc = 0.7. 



De Sitter causet 



The causet in Fig. HI is generated by sprinkling a number of points over a patch in the 1 + 1-dimensional de Sitter 
spacetime, and connecting each pair of points if they lie within the other's light cones. The point density is uniform 
in the de Sitter metric, and the size of the patch is such that the size of the generated causet and its average degree 
are close to those in the considered examples of complex networks. 

In conformal time coordinates, see Section B 1 below, each spacetime point has two coordinates, spatial G [0, 27r] 
and temporal rj g (— 7r/2, 7r/2), with 77 = — 7r/2 and rj — 7r/2 corresponding to the past and future infinities respectively. 
The spacetime patch that we consider is between rj = and rj = tiq > 0, where 770 is determined below. This patch 
is illustrated in Fig. l2ja). To sample N points from this patch with uniform density, we sample N pairs of random 
numbers: N spatial coordinates 9 drawn from the uniform distribution on [0,27r], and N temporal coordinates rj 
drawn from the distribution 



PivlVo) 



sec^ry 
tanryo 



(A2) 



Two spacetime points with coordinates {ri,0) and (r/,0') are then connected in the causet if A9 < A77, where 
A9 — IT — \tt — \9 ~ 9'\\ and A77 = \'>1~ 'n'\ are the spatial and temporal distances between the points in this coordinate 
system. 

To determine 770 we first note that since the point density is uniform, the number of points N is proportioiral to the 
patch volume, and t he pr oportionality coefhcient is a constant point density 5. The volume of the patch is easy to 



calculate, see Section B 1 where we can also calculate the average degree k in the resulting causets, so that we have: 



TV 
k = A5a 



k 

N 



2n6a tan 770, 

2 f Vq 
V tan 7/0 
2 770/ tan 770 + In sec 770 — 1 



In sec 770 — 1 



tan 770 



(A3) 
(A4) 

(A5) 



10 



Q. 

c 
g 

JO 

'i_ 

w 

T3 
CD 

D) 
CD 
T3 

I 

C 



lU 4 




^d = 1 






^d = 3 


10 


'^^. V5yx 




10-^ 
10-^ 


P(k) ~ k"^ ^^^..^^^A 


in-8 


"~ 



10 10 10 10 

node in-degree k. 



10' 



FIG. 6: The in-degree distribution of causets approximated by patches of 1+1 and 3+1 dimensional de Sitter spacetime. The 
causet sizes are A'^ = 23739 and N = 23732, and their average degrees are k = 5.65 and k = 5.25 respectively. 



where a is the spacetime's pseudoradius determining its curvature. Given a target average degree k and number of 
nodes N in a causet, their ratio determines rjo via the last equation. Sampling N points from this patch will then 
yield causets with expected average degree k. The average degree in generated causets will not be exactly equal to 
k since there will be nodes of degree and Poissonian fluctuations of the numbers of nodes lying within light cones 
around their expected values. In Fig. |4j the exact number of sampled points is A^ = 24586 and target value of the 
average degree is A: = 5.53, so that the above equations yield ryo = '^/'^ — 3.86 x 10~^ and Sa^ — 0.151. The resulting 
number of nodes excluding nodes of degree and their average degree in the generated causet are reported in the 
caption of Fig. |4] 

If we direct links in the causet from the future to the past, i.e. from nodes with larger 77 to nodes with smaller 77, 
then the distributions of in-degrees ki and total degrees k — ki + kg are similar. Fig. [5ja), because the distribution 
of out-degrees ko is not particularly interesting, close to Poissonian, Fig. l5jb). The semantics of this link direction in 
causal sets is similar to that in preferential attachment, where it is often convenient to consider links oriented from 
new to old nodes |461 147) . In preferential attachment the out-degree distribution is not particularly interesting either, 
and is given by the delta- function 5{m), where m = k/2 is the number of connections that each new node establishes. 
Figure |4] shows the degree distribution in the undirected causet, i.e. the distribution of total degrees k = ki + ko, 
since all the other networks shown in the figure are undirected (we consider only reciprocal trust relationships in 
the trust network). As a side note, the direction of links in different directed networks may have very different and 
network-specific semantics. In gene regulations |48j or the considered web of trust, for example, link directions show 
what genes regulate what other genes or who trusts whom, respectively, which may have little in common with the 
temporal link direction in preferential attachment or causal sets, showing what nodes are newer or older. 

In higher dimensions, pre-asymptotic effects become more prominent in causets of similar size and average degree. 
In particular, the exponent of the degree distribution is slightly below 2, see Fig. |6] comparing the causets in the d = 1 
and d — 3 cases. The former causet is the same as in Figsj4K45] while the latter is obtained by a procedure similar 



to the one described above, except that instead of Eqs. (A2|A5), we have 
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and each point has two additional angular coordinates, 9i an d 62 , which are random variables between and it drawn 

The spatial distance A9 between pairs of points 



Bl 



from distributions (2/7r) sin^ 6'i and (1/2) sin 6*2, see Section 

is then computed using the spherical law of cosines, and two points are causally linked if A0 < Arj as before. In 

Fig. [ej the d = 3 causet has target N = 25441 and k = 5.29, yielding 770 = 7r/2 - 4.11 x 10"^ and Sa"^ = 0.267. The 
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FIG. 7: Out-degree distribution in the simulated universe. 



maxiiTLum-likclihood fit of the degree sequence in the d = 3 and d = 1 causets yields 7 — 1.65 and 7 = 1.90, while the 
least-square fit of the complementary cumulative distribution function for node degrees yields 7 = 1.77 and 7 = 1.98. 



5. Simulating the universe 



The results in Section [C] below allow us to simulate the universe causet similarly to the de Sitter simulations 
described in the previous section. The main idea behind the universe simulations is that we can use the exact results 
for a flat universe, but apply them to a closed universe with an arbitrary finite number of nodes in its causet, since 
the real universe is almost flat, and since the degree distribution in both cases is the same. 

Specifically, the causet in Fig. [3] is generated as follows. Unlike the previous section where only conformal time is 
used, here it is more convenient to begin with the rescaled time r = t/a. The current measurements of the universe 
To — to/a — (2/3) arcsinh y^fl^/flM = 0.8458 as the best estimate for this rescaled time in the 



C3 



yield, see Section^ w^.^^_^ 

universe today. According to Eqs. (C36|C37), the scale factor is 



R{t) = asinh^ -r. 



(AlO) 



where a is a free parameter that we can set to whatever value we wish, since the degree distribution does not depend 
on it, see Section fC 5[ We wish to set a to the value such that the generated causet would have a desired number of 
nodes 



N ■■ 



1 



n 6a (sinh3ro — 3to) , 



(All) 



where S, the node density, is yet another free parameter that docs not affect the degree distribution according to 
Section C5 The scale factor in Eq. (AlO) means that the temporal coordinate that we assign to each of these N 



nodes is a random number r G [0, tq] drawn from distribution 

smh Stq — 3to 
Having times t assigned, we then map them, for each node, to conformal times rj via 



.23, 
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sinh^ 



(A12) 



(A13) 



The spatial coordinates 0, 9i, and 62 are then assigned exactly as in the previous section, and the causet network 
is also formed exactly the same way, i.e. by future— J-past linking all node pairs whose temporal distance Ary exceeds 
their spatial distance A9. 
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Figure [s] shows the in-degree distribution for a causet generated with a = 2.01, a = 1, and S = 10^. The resulting 
number of nodes in the causet is iV = 10®. The analytic solution curves are obtained using the approximations for the 
in-degree distribution derived in S ection [C| The key equations are Eq. (C21 1, yielding the Laplace approximation for 
the in-degree distribution, and Eq. (C42 1 expressing conformal time as a function of the current value of the scale factor. 



The precise steps to numerically compute the analytic solution for the in-dcgrec distribution are listed in Section C 7 



The perfect match between the simulations and analytic solution in Fig. [3] confirms that the approximations used in 
Section [O to derive the analytic solution yield very accurate results. 

The out-degree distribution in the same simulated causet is shown in Fig. [Tj The analytic solution is obtained by 
numeric evaluations of Eqs.( C13|C17 |. The distribution appears to be uniform over a wide range of degree values, 
which is quite different from the Poissonian out-degree distribution in the sparse de Sitter causet in Fig. [Mb) because 
j^d-i-i _ 2^q4 ^ 1 i]2 the universe simulations, whereas Sa'^'^^ < 1 in the de Sitter simulations in the previous section. 
Depending on whether Sa'^'^^ is smaller or larger than 1, (asymptotically) de Sitter causets have either power-law in- 
degree distributions with j — 2 and Poissonian out-degree distributions, or double power-law in-degree distributions 
with 7 = 3/4 and 7 = 2, and non-trivial out-degree distributions of the type shown in Fig. [t] 
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Appendix B: Asymptotic equivalence between causal sets in de Sitter spacetime and complex networks in 

hyperbolic space 

1. De Sitter spacetime 

The d + 1-dimcnsional dc Sitter spacetime |26| is the exact solution of the Einstein equations 

G^. + A5^, = (Bl) 

for the empty universe with positive vacuum energy density. In these equations, G^^ = i?^^ — i^Rg^i, is the Einstein 
tensor, R^^, the Ricci tensor, R the scalar curvature, g^n the metric tensor, A the cosmological constant, and all the 
notations are in the natural units with the speed of light c = 1. This spacetime can be represented as the one-sheeted 
d + 1-dimensional hyperboloid of constant positive scalar R and Gaussian K curvatures [49] 



- z^ + zf + . . . + z^+i =a'= ' ^^ = ^ = - (B2) 

borrowing its metric from the d + 2-dimensional ambient Minkowski space with metric 



ds'^ = -dzl + dzl + ...+ dzj^^. (B3) 

This hyperboloid in the d = 1 case is visualized as the outer hyperboloid in Fig.|2] As a side note, the curvature of the 
same hyperboloid in the Euclidean metric is everywhere negative but not constant. For d = 3, the Hubble constant iJ, 
vacuum energy density pA, cosmological constant A, and the hyperboloid pseudoradius a, scalar curvature R, and 
Gaussian curvature K are all related by 

where G is the gravitational constant. 

De Sitter spacetime admits different natural coordinate systems with negative, zero, or positive spatial curvatures, 
which are not to be confused with the positive curvature of the whole spacetime. Here we use the standard coordinate 
system (t, 0i, . . . , 6d) with positive spatial curvature that covers the whole spacetime: 

Zf) = a sinh — , (B5) 

a 

zi = a cosh - cos 6*1 , (B6) 

a 

Zd — acosh - sin^i . . .sin0(i_i cos^d, (B7) 

a 

Zd+i — acosh — sin^i . . . sin0rf_i sin^rf, (B8) 

where i e M is the cosmological time of the universe, and 6*1, ... , 9d-i € [0, tt] and Od G [0, 27r] are the standard angular 
coordinates on the unit d-dimensional sphere S'^. The time t at spacetime point P is also the Minkowski length of the 
arc connecting P to the corresponding point at time i = belonging to the zq = slice of the hyperboloid, see Fig. [2] 
In these coordinates the metric takes the Friedmann-Lemaitre-Robertson- Walker (FLRW) form for an exponentially 
expanding homogeneous and isotropic universe with positive spatial curvature: 

ds^ = -dt^ + a^ cosh^ - [dd"^ + sin^ 9 dnl_^) = -di^ + a^ cosh^ - dfl^, where (B9) 

dnl ^ del + s\r?eidel + ... + sm^ei...s\T?ed-idel (bio) 

is the metric on S''. That is, at each time t, the universe is a sphere of exponentially growing radius acosh(t/a) and 
volume 

/ f X"* 2 ^ 

V = (Jd \ a cosh — I , where ad = — ,,,.,, (Bll) 
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FIG. 8: Causal structure of the d + 1-dimensional de Sitter spacetime. d — 1 dimensions are suppressed, so that each point 
represents ad— 1-sphere. AS is the distance between A and B on the d-dimensional unit sphere S"*. 

is the volume of S''. Figure [2] visualizes this foliation for d = 1, in which case these time-slice spheres are circles, and 
their volume is the circle's circumference. 

To study the causal structure of de Sitter spacetime, it is convenient to introduce conformal time 77 e (— ^, ^) via 



sec?7 = cosh — . 
a 

These conformal-time coordinates are convenient because the metric becomes 

ds^ = a^ sec^ 77 {-drf- + dd^ + sin^ d d^\_^ , 



(B12) 



(B13) 



so that the light cone boundaries defined by As — are straight lines at 45° with the (77, 6*) axes, see Fig.lsj Therefore, 
in this figure, point A at time 77 S [0, r]') lies in the past light cone of point B at time 77' G (0, 7r/2) if the angular 
distance A0 between A and B on S'' is less than the conformal time difference r( — r\ between them: 



/ ft. 

A0 < 7/ — 77 = arcsec cosh arcsec cosh — , or approximately 

a a 



A0 < 2 e 



(-^-^) 



2e-^, 



(B14) 

(B15) 



where the last approximation holds for t' ^ i ^ 1, and where we have used the approximation rj = arcsec cosh(i/a) 
7r/2 — 26"*/", which is valid for large t. 

The volume form on de Sitter spacetime in the conformal and cosmological time coordinates is given by 



dV = (a sec 77) drjd^d— (a cosh- | dtd^, 



"""(l) 



c^dtd^d, where 



d$« 



sin'^"^ 01 sin'*"^ 62... &in6d-idei d92 . . . dOa 



(B16) 
(B17) 



is the volume form on S'', i.e. / d^d — Cd. 



2. Hyperbolic space 



The hyperboloid model of the d + 1-dimensional hyperbolic space 1361 is represented by one sheet of the two-sheeted 
d -\- 1-dimensional hyperboloid of constant negative scalar R and Gaussian K curvatures 
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(B18) 



borrowing its metric from the d -\- 2-dimensional ambient Minkowski space with metric ( B3 ) . This hyperboloid in the 
d = 1 case is visualized as the inner hyperboloid in Fig. [2] As a side note, the curvature of the same hyperboloid in 
the Euclidean metric is everywhere positive but not constant. The standard coordinate system (r, 6'i, . . . , 9d) on this 
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hyperboloid is given by 

zo = 5 cosh ^, (B19) 

zi = 5sinh^cos6li, (B20) 



Zd 
Zd+l 



T 

= 6sinh - sin6'i . . . sin6'£(_i cos6'(i, (B21) 



= 5sinh - sin^i . . . sin0(i_i sin^d, (B22) 

where r E M+ is the radial coordinate, and 6i, . . . , Od^i E [0, n] and 9d E [0, 27r] are the standard angular coordinates 
on the unit d-dimensional sphere S . The radial coordinate r of point P on the hyperboloid is the Minkowski length 
of the arc connecting P to the hyperboloid vertex, which is the bottom of the inner hyperboloid in Fig. [2] As a 
side note, the Euclidean length te of the same arc is given by the incomplete elliptic integral of the second kind 

rE = —iE{ir, 2). 

3. Hyperbolic model of complex networks 

In the hyperbolic model of complex networks .35i , networks grow over the d + 1 dimensional hyperbolic space of 
Gaussian curvature K = —1/6^ according to the following rule in the simplest case. New nodes n are born one at a 
time, n = 1,2,3,..., so that n can be called a network time. Each new node is located at a random position on S''. 
That is, the angular coordinates (0i, . . . , 6d) for new nodes are drawn from the uniform distribution on S'*. The radial 
coordinate of the new node is 

h Ti 

r = 2-ln-, (B23) 

a V 

where v is a, parameter controlling the average degree in the network. Upon its birth, each new node connects to all 
the nodes lying within hyperbolic distance r from itself. In other words, the connectivity perimeter of new node n 
at time n is the hyperbolic ball of radius r centered at node n. The hyperbolic distance x between two points with 
radial coordinates r and r' located at angular distance A0 is given by the hyperbolic law of cosines |36j : 

f r r' r r' \ A0 

X = b arccosh cosh - cosh -, sinh y sinh — cos A9 « r + r' + 26 In — . (B24) 

\ b b b b J 2 

Therefore new node n' connects to existing nodes n < n' whose coordinates satisfy 

X < r' , or approximately (B25) 

An 

r + 2bln— < 0. (B26) 

This construction yields growing networks whose distribution P{k) of node degrees fc is a power law, P{k) ^ k~'^ , 



with 7 = 2. Indeed, according to Eq. (B23), the radial density of nodes at any given time scales with r as p{r) ~ e"'', 
where a — d/{2b). One can also calculate, see ^50], the average degree of nodes at radial coordinate r, which is 
k{r) ~ e"'^'", where f3 = a. The probability that a node at r has degree k is given by the Poisson distribution with 
the mean equal to k{r). Taken altogether, these observations prove that 7 = a//3 + 1 = 2. The networks in the model 
also have strongest possible clustering, i.e. the largest possible number of triangular subgraphs, for graphs with this 
degree distribution, and their average degree is 

k « 2'^+^—iylnn, where v^ = ^^^ (B27) 

ad a 

is the volume of the unit d-dimensional ball. The model and its extensions describe the large-scale structure and growth 
dynamics of different real networks, e.g. the Internet, metabolic networks, and social networks, with a remarkable 
accuracy [55] . 
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4. Duality between de Sitter causal sets and complex networks 

To demonstrate the asymptotic equivalence between growing networks from the previous section and causal sets 
growing in de Sitter spacetime, we: 

1. find a mapping of points in the hyperbolic space 13''+-'^ to de Sitter spacetime dS'^'^^ such that: 

2. the hyperbolic ball of new node n e H''+^ is asymptotically identical to its past light cone in (i§'^+^ upon the 
mapping, and 

3. the distribution of nodes after mapping is uniform in dS'^'^^. 

a. Mapping 

A mapping that, as we prove below, satisfies the properties above is remarkably simple: 

t = r, (B28) 

a ^ 2b. (B29) 

That is, we identify radial coordinate r in M'^+^ with time t in rfS''+^, keeping all the angular coordinates the same — see 
Fig. [2] for illustration. 

b. Hyperbolic balls versus past light cones 

The proof that the hyperbolic balls of new node connections map asymptotically to past light cones is trivial: 
inequalities (|B26) and (B15) are identical with the mapping above. Figure ^d-f) visualizes the approximation 



accuracy at different times. 



c. Uniform node density 



The proof that the node density after mapping is uniform in dS'^'^^ is trivial as well. Indeed, according to Eq. (B23) 
with 26 = a, network time n is related to cosmological time t via 

n = i/e»*, so that (B30) 

dn = v-ei^dt. (B31) 

a 



Using the last equation, we rewrite the de Sitter volume element in Eq. (B16) as 



26''+i 

dV ^ —dnd^^d- (B32) 

vd 



By construction in Section B 3 the node density on S'' is uniform and equal to 1. Therefore the number of nodes dN 
in clement dn d^d is 



Combining the last two equations, we obtain 



dN=—dnd<^d- (B33) 



dN = 5dV, where (B34) 

"^ (B35) 



2¥+^ad 
meaning that nodes are distributed uniformly in d^'^^^ with constant density 5. thus completing the proof. 
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A causal set growing in de Sitter s pace time with node density 6 in Ecl( B35 ) is thus asymptoticaUy equivalent to a 
growing complex network in Section B3 with average degree k in Eq. (B27|. Combining these two equations, we can 
relate 6 and k to each other: 



k « 2vd5a'^t. 



(B36) 



An important consequence of this asymptotic equivalence is that the degree distribution in both cases is the same 
power law with exponent 7 = 2. 



Lorentz invariance 



Ijd+l 



IS 



Here we show that the described duality is Lorentz- invariant. The group of isometrics of dW^'^^ and 
SO{d + 1, 1), i.e. the Lorentz group of the ambient Minkowski space M''+^. Any g G SO{d + 1, 1) preserves volumes, 
maps light cones in dS'^'^^ to light cones, and sends balls in H^+^ to balls of the same radii. The fact that the mapping 



$ : d§'*^ 



i-> 



ird+l 



described in the previous section is invariant under the SO{d + 1) subgroup of SO{d + 1, 1) is 



e SO(d + 1, 1). This element defines a new coordinate system in 
j, where gij is the matrix representation of g, e.g. if g is the Lorentz boost with rapidity (j) in the 



immediately obvious. Consider now any element g 

— — sinh 



T hese new z-coord inates 
|B8|b19|B22| . The 



zi direction, then goo = 5ii — cosh 0, gi2 — g^i — — sinh 0, and gkk = 1, 5fci = for fc, ^ > 1 
define new coordinates (i, ^i, . . . , 9d) and (f , ^i, . . . , 0d) on dS''''"^ and 11'*+^ via the same Eqs. (B5 
same g, considered as an active transformation, sends the patch between t = and t = t' in dS +^ in Fig.l2Jto another 
patch. As a side note, since SO{d+ 1, 1) is not compact, there exists no compact subset of (i§'^'+^ (or M''^) invariant 
under the action of SO{d + 1, 1). Yet one can check that g sends the patch between t — {) and t = t' to the patch 
between t — Q and t = t', the past light cone of point P aX t = t' — to the past light cone of P's image P ai t — t' , the 
ball of radius r' centered at r = in EI''+^ — to the ball of the same radius centered at f = 0, a nd so on. Therefore, 
mapping 4> : dS'^'^^ i— >■ EI''+^ after the transformation induced by g is exactly the same as in Eq. (B28I, i.e. t — r, and 
under the action of SO{d + 1, 1) this mapping transforms as $ = g o $ o g~^. 
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Appendix C: The universe as a causal set 

De Sitter spacetime is the spacetime of a universe with positive vacuum (dark) energy density and no matter or 
radiation. Since the real universe does contain matter, its spacetime deviates from the pure de Sitter spacetime. At 
early times, matter dominates, leading to the Big Bang singularity at i = that de Sitter spacetime lacks. At later 
times, the matter density decreases, while the dark energy density stays constant, so that it starts dominating, and the 
universe becomes asymptotically de Sitter. The universe today is at the crossover between the matter-dominated and 
dark-energy-dominated eras, since the matter and dark energy densities pm and pA are of the same order of magnitude, 
leading to rescaled cosmological time r — t/a ^ 1 — the so-called "w/ij/ now?^ puzzle in cosmology pSHST] . 

To quantify these deviations from pure de Sitter spacetime, and their effect on the structure of the causal set of the 
real universe, we calculate its degree distribution in this section. This task is quite challenging, and in what follows 
we first provide the exact analytic expression for the degree distribution, and then derive its approximations based on 
the measured properties of the universe. These approximations turn out to be remarkably accurate because according 
to the current measurements, the universe is almost flat. 

We emphasize that the approximations in this section are based on the exact solution of the Einstein equations for 
a flat universe containing only matter and constant positive vacuum energy, i.e. constant cosmological constant A, 
and that we assume that the universe is governed by this solution at all times. This assumption is a simplification of 
reality for a number of reasons. For example, there are a plenty of cosmological scenarios, such as eternal inflation j51j . 
in which the ultimate fate of the universe deviates from the asymptotic solution that we consider. There are also a 
variety of other models with non-constant A, yet the standard ACDM model with constant A is a baseline describing 
accurately many observed properties of the real universe |52| . partly justifying our simplifying assumption. We also 
note that by relying on the exact solution for a universe containing only matter and A, we effectively neglect the 
earliest stages of universe evolution such as the radiation-dominated era or inflationary epoch. There is no consensus 
on how exactly the universe evolved at those earliest times. Since the radiation-dominated era ended soon after the 
Big Bang [S3], these details are unlikely to have a profound effect on the universe causet's structure at much later 
times. 

1. Exact expression for the degree distribution 

The metric in a homogeneous and isotropic universe takes the FLRW form: 

1 



ds^ = ^dt^ + R\t)idx + —sm:' y/Kx (sin^ 9 dO^ + dcj)^) \ , wheve (CI) 

{sin X \i K = 1 (closed universe with positive spatial curvature), 
X^ \i K ^ (flat universe with zero spatial curvature) , (C2) 

sinh X if ^ = — 1 (open universe with negative spatial curvature). 

In the above expression, coordinates e [0, 27r] and 9 G [0,7r] are the standard angular coordinates on S^, while 
X S [0, 7r] if X = 1, or x G [Oi oo\\i K — 0, —1, is the radial coordinate in a spherical, flat, or hyperbolic space. Finally, 
t £ [0, oo] is the cosmological time, and R{t) is the scale factor, finding an appropriate approximation to which in 
the real universe is an important part of our approximations in subsequent sections. In this and the next sections, 
however, all the expressions are valid for any scale factor. 

As with pure de Sitter, it is convenient to introduce conformal time 77, related to cosmological time t via 

In conformal time coordinates, the metric and volume form become 

ds^ = R^{ri)[-dri^ + dx^ + —sm^ VKx {sin^ 9 dO"^ + dcj)"^)^ , (C4) 

dV = —R'^{ri)s\u^\\fKx\^'^^Sdridxded(t>. (C5) 

K L J 

As with pure de Sitter, we orient edges in the universe's causet from future to the past. Therefore if 77 is the current 
conformal time, then the average in- and out-degrees fci_o(^'|^) of nodes born at conformal time 77' < 77 are simply 
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proportional to the volumes of their future and past light cones Vf,p{ri'\r]): 



(C6) 

(C7) 



where the coefhcient of proportionality is the Planck-scale node density in the spacetime, which we take to be the 
inverse of the Planck 4-volunic: 



1 



1.184 X 10' '^ s"*, where 



tp = 5.391 X 10"^"* s 

is the Planck time. 

Thanks to conformal time coordinates, the expressions for volumes Vf^p{r]'\ri) are easy to write down: 



(C8) 
(C9) 






Kx 
Kx 



de &in6 

JQ 

dO sin 

"'0 



27r 



Computing the three inner integrals, we obtain 



ktiil'lv) = ^ 



koiri'lv) = ^ 



K 



K 



1-1 ' 



2(V' - V') 



2(r?' - v") 



K 

1 



sm 



K 



2^/K{v"-^^R\v")drJ", 
2\/i^(?7'-r;")lli?*(ry")d?y". 



(CIO) 
(Cll) 

(C12) 
(C13) 



As shown in [20' , Lorentz invariance implies that nodes/events of the causet are distributed in spacetime according to 
a Poisson point process. Therefore, as explained in |54j . to find the in- or out-degree distributions P(k,ri) at time 77, 
we have to average the Poisson distribution 



which is the probability that a node born at time 77' has degree k, with the density of nodes born at time 77' 

U{v) = rR\v')dv' 
Jo 



is the time-dependent normalization factor. The result is 

P(fc,77) = / p{k\r)',r])p{r]'\r])dr]' = _^yY^^ / 



'k{v'\v)f e-'''^'^'^'^^ R\v') dr,' . 



(C14) 

(C15) 
(C16) 

(C17) 



The above expressions are valid for both in-degree (set k = ki and k(r]'\ri) = ki{r]'\ri)) and out-degree (fc = ko, 
k{ri'\rj) = ko{ri'\ri)) distributions. In what follows we focus on the in-degree distribution. 



2. First approximation using the Laplace method 



Equations ( C17|C12 ) give the exact solution for the in-degree distribution with an arbitrary scale factor R{ri), but 
it is difhcult to extract any useful information from these expressions, even if we know the exact form of R{ri). The 
first step to get a better insight into the causet properties is to rewrite ki{ri'\T]) as 



kiiv'lv) — ki{Q\v)P{v'\v)^ where F{ri'\f]) 






(CIS) 
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The in-degree of the oldest node ki{0\T]) is an astronomically large number, both because the node is old, so that 
its future light cone comprises a macroscopic portion of the 4- volume of the whole universe, and because the degree 
is proportional to huge S. However, function F{r]'\r]) is a monotonically decreasing function whose values lie in the 
interval [0,1]. We can thus use the Laplace method to approximate the integral in Eq. (C17) in the limit S ^ 1. 
Introducing •q*{ki), which is the solution of the transcendent equation 



and function 



$(fc„?7) 



the result of this Laplace approximation reads: 



h{v*\v), 



dki{r]'\T]) 



dr\' 



(C19) 



(C20) 



ri'=ri*(ki) 



P{h,v) = < 



5\ R^{tj) 



1 V^^zk^e-^^ R''[v*{h)] 



1 RHv*ik^)] 



l^f{v) 



k,\ 



<i>ih,r]) Af{r]) $(fc,,?7) 



if fc, = 0, 



if l<fc, <fc,(0|r7). 



(C21) 



where we have also used Stirling's approximation fc! « ^/2nk{k/e)''. We see that the shape of the degree distribution 
is almost fully determined by r]*{ki), and that the dist ribut ion is a fast decaying function for degrees above ki{0\ri). 

The expression for the in-degree distribution in Eq. (C21) is now more tractable and gets ready to accept the scale 
factor R{rj) of the universe. Unfortunately, the exact expressions for the scale factor of a closed or open universe 
with matter and dark energy, although known [5S], re sist an alytic treatment, and so does the integral for the average 



degree kilrj'lr]) in Eq. (C12) that we are to use in Eq. (C21 ). Therefore we next develop a series of approximations to 



the scale factor and average degree, based on the measured properties of the universe. 



3. Measured properties of the universe 



The current measurements of the universe 



that are relevant to us here include: 



flA = :-^ PA e [0.709, 0.741], (dark energy density) 



^M — 



SttG 



Pm e [0.2582,0.2914], (matter density) 



K 






€[-0.0133,0.0084], (curvature density), where 



Hn = — e [68.8,71.6] °^ = [2.23,2.32] x 10"^* s"^ is the Hubble constant, 

i?o s • Mpc 

Rq = R(to) is the scale factor at present time to, i.e. the age of the universe, 

io e [13.65, 13.87] Gyr = [4.308,4.377] x 10^^ s, so that 

A = SH^i^A e [1.06, 1.20] X 10^^^ s^^, is the cosmological constant, and 

^^ 1 



'TT^Hq 



e [5.00, 5.32] X 10 s, is the de Sitter pseudoradius. 



(C22) 
(C23) 
(C24) 

(C25) 

(C26) 

(C27) 
(C28) 

(C29) 



The range of values of ilA, ^m, ^Oj and to are given by the 95% confidence bounds in the universe measurements taken 
from the last column of Table 1 in [52], while the values of Hk come from Table 2 there: third row, last column. The 
matter density J^m is the sum of two contributions: the observable (baryon) matter density (^1^ G [0.0442,0.0474]), 
and dark matter density {ftc € [0.214,0.244]). The radiation density is negligible. The universe thus consists of 
dark energy (« 73%), dark matter (« 23%), and observable matter (« 4%). As a side note, the sum of all fi's 
^ rj e [0.95,1.04] as expected, since ^ fi = 1 is the Einstein/Friedmann equation for a FLRW universe, which we 
recall in the next section. Much more important for us here is that the universe is almost fiat, fix ~ 0. 
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4. Approximations to the scale factor and average degree 

In this section we utilize the approximate flatness of the universe to derive and quantify approximations to the scale 
factor of the universe and average degree in Eq. (C12 ). The main idea is that since ^k is small, we can use the exact 



solution for the scale factor in a flat universe {Qk = 0) with matter and dark energy (positive cosmological constant), 
which is actually quite simple |26) . 

We first recall that the 00-component of the Einstein equations in the FLRW metric is the first Friedmann equation: 



R^ + K 

i?2 



8 



vrG 



PA 



This equation can be rewritten ,56j in the following form: 

cR/Ro 



Ro 
~R 



dx 



Phi 



(C30) 



Assuming non- vanishing positive dark energy and matter densities (J^a > and i7j\/ > 0), using rescaled time 

t 



a 



and introducing two parameters 



Ro\Iyv~' ^'^"^ 



n 



K 



y^A^M 



e [-0.0368,0.0232], 



the integral in Eq. (C31) simplifies to 



dx 
/o x\n-\-ex^ + x^ 

Solving it for r = r{T,e), we conclude that the solution for the scale factor becomes 

R{t) = ar(r;e). 



(C31) 

(C32) 

(C33) 
(C34) 

(C35) 
(C36) 



This is definitely not the only way to write down the scale factor in terms of the parameters of the universe. Yet 
written in this form, the rescaled scale factor r(r; e) is a dimensionless function of its dimensionless arguments, and 
parameter e is close to zero, so that ^(r; e) « r(T;0), where r(T, 0), i.e. the rescaled scale factor for a flat universe 
with positive matter and dark energy densities, is quite simple: 



r{T; 0) = sinh^ -r. 



(C37) 



Figure ^left) shows a good agreement between the numerical solution for r{T; e) and the analytical expression for 
r{T; 0) in Eq. ( |C37 ) for the range of the values of e allowed by measurements in Eq (C34). Therefore in what follows 
we can use Eq. (|C37), even if the universe is "slight ly clo sed" or "slightly open." 



To approximate the average degree ki(ri'\ri) in Eq. ( C12 ), we first compute, using the scaling Eq. ( C36 ), the conformal 
time 



V = 



dt' 
W) 



dr' 



r{T';e) 



r(r;e) 



dx 



'VT^ 



(C38) 



at the future infinity t — t — r(T]e) = oo. For the extreme values of e allowed by measurements, we obtain 
^700 = 0.5403 for the closed universe with e ~ —0.0368 and rj^ = 0.4260 for the open universe with e = 0.0232. If 
e — > 0, then we also have ?7oo — ^ 0, yet if e is exactly zero, then rjoo is an undefined constant, since in the fiat case, Rq 
is a free parameter that can be set to an arbitrary value without affecting anything. If the universe is not exactly fiat. 
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FIG. 9: Left: Relative error between functions r(T\e) and r(T;0), |r(r;e) — r(r; 0)|/r(r; e), for the two extremes of the 95% 
confidence interval of parameter e in Eq. (C34l, as a function of rescaled time r. In the worst case, the relative error is 



around 1%. Right; Relative error of approximation in Eq. (C39l for x £ [0, 7700] 



then the above values of 1)00 allow ed by measurements are well below 7r/2, so that we can safely replace expression 



Ix — sm{2\fKx) I \rK jK in Eq. (C12) by its Taylor expansion around zero, 

1 



1 



2x 



s\TL{2y/Kx) 



(C39) 



Indeed the maximum possible relative error between the left and right hand sides in the last equation &i x — rjoo is 
around 6%, see Fig. [91 right), whereas such error is zero for an exactly flat universe with e = K = 0, since the left and 
right hand sides are equal in this case. Therefore the approximation to the average degree in Eq. (C12) 



\^ jyii 



Hl'\v) = \^5 j\^"~^rR\i')dv" 



(C40) 



is exact for flat universes, and almost exact for slightly closed or open universes with e within the range of Eq. (C34| 



5. Scaling relations 



To proceed to our final approximations to the degree distribution in the next section, we derive some scaling relations 
that tremendously simplify the calculations, provide important insights into the properties of the degree distribution, 
and suggest methods to simulate the causal set of the universe. 

We begin with the scaling relation for conformal time 77. Let us introduce the rescaled conformal time 



a 



V 



(C41) 



and let notations y ~ fi.2....i^) mean "y is a function of x." Given the definition of conformal time in Eq. (C38), this 
rescaled conformal time is a function of rescaled cosmological time r: 



C = /i(r) 



r{T]e) 



dx 



'-VT 



-3 e-+0 



>2v/K7yO)2Fi(^,i ^,-r3(r;0)). 



6' 2' 6' 



where r(T;0) is given by Eq. (C37), and 2F1 is the hypergeometric function. This scaling means that 

which in turn implies that the average in-degree of nodes born at time 77' 

^il??'!??) = 9/3 (C', C) = qfi (t', t) , where 
q = 5a , 



(C42) 

(C43) 

(C44) 
(C45) 
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while the average in- or out-degree in the causet 

Hv) = Uv) = r W\v)pW\v) dri = qh (C) = 9/6 (r) . (C46) 



Plugging scaling Eqs. (C44|C45| into the expression for the in-degree distribution in Eq. (C21 ), we conclude that this 



distribution scales for fc^ > 1 as 

P{k,,i) = -Q(K,,r), where (C47) 

At. = -, (C48) 

and the rescaled in-degree distribution Q{Ki,T) is some function of rescaled degree k^ and time r. This result is 
important for several reasons: 

1. It defines a characteristic degree q = Sa"^ and characteristic time a such that the degree distribution depends 
only on the rescaled dimensionless variables Ki = ki/q and r — t/a. Therefore we are free to set a — 1 and 



(5 = 1, and the remaining task is to find an explicit form of the rescaled distribution Q{ni, r) in Eq. (C47) 



2. Even though the scale factor in Eq. (C36) does depend on a, the in-degree distribution P{ki, t) does not depend 
on a. This is important from the practical point of view, because we can carry out all the calculations setting 
a = 1 as well. 

3. Having set a = a = i5 = 1, the in-degree distribution explicitly does not depend on any physical parameters of 
the universe. The causets in flat or open universes are graphs with obviously infinite numbers of nodes from the 
very beginning of the universe, while the causets of closed universes are always finite, and the number of nodes 
in them depend on time and on all the three parameters a, a, and 5. Yet the in-degree distributions in all these 
causets at a given time are all the same, and in particular, do not depend on the number of nodes at all. This 
somewhat unexpected result deserves further explanation. We note that parameter a = i?o \/^mI^a depends 



on the physical properties of the universe. In particular, if the universe is closed, it defines, via Rq in Eq. ( C24 1 , 
the radius of the three-sphere, i.e. of the spatial part of the FLRW metric. Since the degree distribution does 
not depend on a, we can freely take the limit a — ^ cx) corresponding to a flat space, without affecting the degree 



distribution. In other words, we are free to choose a closed universe with K — +1 in Eq. (CI), and this choice 
can be considered as a degree-distribution-preserving compactiflcation of an infinite flat or hyperbolic space 
with an infinite number of nodes in it, into a compact spherical space with a finite number of nodes. 

4. The previous point informs us how to simulate the universe causet on a computer. We cannot simulate infinite 
causets for fiat or open universes, but since the degree distribution is the same in closed universes, we are free 
to simulate the latter. One can check that in a closed universe, the number of nodes in a causet grows with 
time as 

N{t) = Saa^fr (r) . (C49) 

We can set a = 1, and since the in-degree distribution does not depend on the number of nodes N, then for 
any given time t = ar = t and for any given value of 5, we can generate graphs with any desired size N by 
choosing the value of a accordingly. The degree distribution obtained from this simulation is then to be re-scaled 



according to Eq. (C47) to obtain the rescaled distribution Q{Ki,T). Keeping fixed t = t and N, we can also 



explore different regions of Q(Ki, r) by changing the value of S. 
We have to emphasize that the points above, and the scaling relations in this section, are valid for non-flat universes 



only when t/oo is well below n/2, so that the approximation in Eq. ( C39 ) is valid for all x £ [0, t/oo] ■ When this condition 
does not hold, e.g. when e is large and the universe is far from being flat, then all these scaling relations break down. 
In particular, we cannot claim that in that case the in-degree distribution would still be the same. The latest 
measurements of the universe parameters indicate that the scaling laws derived in this section do apply to the real 
universe, but even if this were not the case, these laws would still be valid for small values of conformal time. 

6. Final approximations to the degree distribution in the matter- and dark-energy-dominated eras 

We now have all the material needed to derive the final approximations to the degree distribution. According to 



Section C 4 the exact solution for the scale factor in a flat universe is a good approximation for the scale factor in the 
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real universe within the measurement-allowed range of parameter values, so that we will use this approximation here, 
i.e. we set e = 0. Using the scaling results from Section [C 5[ we also set a = a = S = 1, which is equivalent to working 
with rescaled cosmological time r, rescaled conformal time C, rescaled scale factor r{T) = sinh^'"^ (3t/2), and rescaled 
in-degree k^. Since everything is rescaled in this section, we omit word "rescaled" in front of the variable names. 

We consider two eras of the universe evolution. The first era characterizes the causal set of the universe at early 
times, when matter dominates. The second era deals with the aged universe at large times, when dark energy 
dominates. In the limit of infinite time (future infinity), we derive the exact solution for the degree distribution in a 
flat universe. 



a. Matter-dominated era 



In this era with r ^ 1, the scale factor and conformal time can be approximated by 

r{r) = (^rj 

( = 2 (-t) . so that 



r-(C) 



(C50) 
(C51) 
(C52) 



Using these expressions in the rescaled version of Eq. (C40), we obtain the average in-degree at time C of nodes born 
at time C': 



s,(c'|C) = ^c'' 



1 3 C' 3 fCV 1 fC^ ^ 



12 11 C 10 VC/ 9 VC/ 1980 VC 



1 fC 



The maximum average in-degree, i.e. the average in-degree of the oldest nodes, is 

TT 



«.(o|C) 



9216 



C'" = 97rT^ 



(C53) 



(C54) 



We note that if we reinsert constants a and S, and cosmological time t using Eq. (C44), the average in-degree of the 
oldest nodes becomes 



h{0\t) =97r 



(C55) 



where tp is the Planck time, see Eq. (C8|. As expected, this in-degree does not depend on a and, consequently, on 
the cosmological constant A. 



For degrees well below the maximum degree, Ki ^ Qttt^, the Taylor expansion of Eq. (C53l around C' — C yields 

4 



s,(C'IC) 



>12 



3072^ V C 



c 



Using this expression in the rescaled version of Eq. (C21), we obtain 

Q(k*,t) = ^(-) -^ ifK, <9^T^ 

and, after reinserting all the physical constants. 



4 V TT 



pfeo^^m'f.,-'" ifi<..«9,(± 



(C56) 



(C57) 



(C58) 



with a soft cut-off at 9tt {t/tp) . We thus observe that the degree distribution does not depend on the cosmological 
constant A, and that it is a power law P{ki,T) ^ k~'^ with exponent 7 = 3/4. 
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b. Dark-energy-dominated era 



We now analyze the future fate of the universe at r 3> 1, which is shghtly more intricate. To begin, we derive 
a co uple of expressions that we use in simulations and that are valid for any r. We first recall that, according to 
Eq. ( C42 ) , conformal time is related to the current value of the scale factor r via 



C(r) = 2V^ 2F1 



1 1 7 
6'2'6' 



(C59) 



which is a monotonously increasing function of r. Therefore we can use scale factor value r as a measure of time 
instead of (. Using this observation in Eq. (C40), we write the average in-degree at time r of nodes born at time r' as 



Ki{r'\r) = 



47r 



[C(^)-C(r')] 



VT 



r-3 



Similarly, function ^{ni,r) in Eq. (C20) becomes 



$(Ki,r) =47r / 

J r* 



r*(K,) Vl- 



;:^' 



(C60) 



(C61) 



where r*(Ki) is the solution of equation Ki{r*\r) = Ki 



Assuming now that time is large, we see from Eq. ( C60 1 that the maximum average in-degree scales at r 3> 1 as 

dx 47r 



K.(0|r) 

Si 00 



47r 

T 



^'«7f 



= c(cx3) = r 

Jo 



- X 
dx 



-3 



^VT 



.-3 



So 



y^ 



^dx = y {C.cK>rf , 



where 



(C62) 
(C63) 



is the conformal time at the future infinity. Since r grows exponentially with time in this era, so does the average 
in-degree of the oldest nodes according to Eq. (C62). In the long time limit, and for degrees well below the maximum 
average degree, Ki <C Ki(0|r), we have r > r' ^ I. Keeping the first two terms in the Taylor expansion of Eq. (C59 ) 
for r ^ 1, we approximate conformal time by C(r) ~ Coo — ^/f- Inserting this approximation into Eqs. ( C60|C6T |, 
and neglecting x~^ there, we obtain the implicit expression for the in-degree distribution at the future infinity, 



Q(Ki,00) 



9 



1 



4tt x{k,)^[x{k,) - If 



where x{Ki) is the solution of equation 



27r 



[{x - l)(2a;2 -7x + 11) - 61na;] 



(C64) 



(C65) 



in the region x > 1. The last two equations give a nearly exact solution for the asymptotic in-degree distribution in 
a flat universe, because all the approximations that we have made so far become exact in the i,r — > 00 limit. The 
only approximation that is not rigorously exact is the one due to the Laplace method in Section [C 2[ yet given the 
astronomical value of node density S in Eq. ( |C8[ ) used in this approximation, it can be also considered exact. We also 
note that the distribution is properly normalized since J. Q{x,oo) (kOk dx = 1. 

Fi nally , we find approximations to this exact solution for small and large degrees. If Ki <^ 1, then the solution of 
Eq. (C65) scales as x{K i) = 1 + (3Ki/7r)^/'*, whereas for Ki :$> 1, the scaling is x{Ki) = (9Ki/47r)^/^. Substituting these 
scalings into Eq. (C64), and neglecting insignificant terms there, we obtain 
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(C66) 



where the crossover degree value c = 1.66 is given by equating the two approximations above. Figure [TO] shows the 
exact numerical solution for the in-degree distribution at t = 00, and juxtaposes it against these two approximations. 
The match is remarkable. 
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FIG. 10: Exact numerical solution versus approximations for the in-degree distribution in the universe causet at the future 
infinity. The solid line shows the numeric solution in Eqs. ( C64|C65 l, while the dashed (-dotted) lines are the two approximations 
in Eq. ( C66 1 for Hi <^ c and Ki ^ c. The crossover degree value c is shown by the dotted vertical line. 



7. Numeric evaluation of the analytic solution 

Given the results in the previous sections, the precise steps to numerically evaluate the analytic solution for the 
in-degree distribution in the universe at any given time are: 

1. For a given rescaled age of the universe r, e.g. r = tq, compute the current value of the rescaled scale factor 



r^r) — sinh^ -r, 



(C67) 



and the normalization factor 



A/'W = TT(sinh3r-3r). 



(C68) 



2. According to Eq. (C42), define function 



C(0 = 2Vf.^(| 1,1,^.3 



(C69) 



where 2-^1 is the hypergeometric function. 



3. Generate a (log-spaced) sequence of rescaled in-degrees k^, and for each value of Hi in the sequence, find 
numerically the solution r*{Ki) G [0, rir)] of equation 
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4. With r*{Ki) and r(T) at hand, compute numerically the integral 

r^'^^ T x^dx 

^K,, r{T)) = in / [C(x) - C{r*{n,)f -—^ 



(C71) 



5. Finally, according to Eq. (C21 ), the value of the rescaled in-degree distribution at ni is given by 
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8. Universality of 7 = 3/4 scaling for small degrees 

The portion of the degree distribution with exponent 7 — 3/4 for small degrees below ~ Sa'^ remaining in the 
universe causet even at long times may appear as a paradox at the first glance. Indeed, at long times, the universe is 
in its accelerating era dominated by dark energy, and the scale factor grows exponentially, versus polynomial growth 
in the matter-dominated era. Yet the degree distribution for small degrees behaves exactly the same way as in a 
matter-dominated universe, i.e. it has the same exponent 7 = 3/4. 

The intuitive explanation of this paradox is as follows [52| . The degree distribution in the range of small degrees 
is shaped by spacetime quanta born at times 77' near the current time rj. The future horizon radii of these nodes are 
smaller than the Hubble radius I/Hq ^ a. Therefore these nodes do not yet "feel" the acceleration of the universe. 
For them the universe expands as if it was matter-dominated. 

To formalize this intuition we check analytically that exponent 7 = 3/4 for degrees fc^ <C Sa^ is universal for any 
scale factor. Keeping only the first term in the Taylor expansion of ki(ri'\ri) at ry' — rj, we approximate 

h{v'\v)^l7rSR*iv)iv'~v)'- (C73) 



Using this approximation to solve equation ki = ki{ri'\r]) for ri*{ki), and substituting the solution into Eq. (C21), we 
obtain 

The degree distribution for small degrees is thus a power law with universal exponent 7 — 3/4, while the scale factor 
is relegated to normalization. The scale factor is important only for old spacetime quanta with large degrees, where 
its exponential growth in asymptotically de Sitter spacetimes is responsible for the emergence of exponent 7 = 2. 
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Appendix D: Related work 

The idea of replacing continuum spacetimc with a graph or network appears in many approaches to quantum 
gravity, as it is a natural way to describe a discrete geometry. Causal Dynamical Triangulations (CDTs), for example, 
are formulated in terms of a simplicial triangulation of spacetime, which can be regarded as a (fixed valence) graph 
in which adjacent simplices are connected by an edge |58H60) . Another popular approach to quantum gravity is Loop 
Quantum Gravity, where the quantum states of geometry are described naturally in terms of spin networks, which 
are graphs embedded into a three dimensional manifold. The edges and vertices in these graphs are colored with 
various mathematical structures, see |61j and references therein. There are a multitude of other approaches whose 
mathematical formulations have a similar discrete network-like character, such as Wolfram's evolving networks |62j . 
"quantum graphity" [63], D'Ariano's causal networks [64j, Requardt's lumpy networks [65], etc. 

Most descriptions of quantum spacetime geometry in terms of a graph structure are of a purely spatial character, 
in that the description of spacetime arises from temporal evolution of the graph whose edges usually connect spatially 
nearest neighbors. The intuition behind this idea is clear, being similar to the idea of approximating a continuous 
space by a fine lattice embedded in it. However, the precise physical meaning of the "time" in which the network 
evolves, and the manner in which the Lorentz invariant nature of special relativity can emerge from such discreteness, 
are often unclear [20] • Causal sets and CDTs are different in that they are described in terms of a graph-like structure 
which has a fundamentally spaceiijTie character. 

The degree of a node is an effective measure of its age in growing complex networks. This simple observation led us 
here to establishing a connection between these networks and the causal networks of discrete spacetime. The notion 
of degree as a measure of time is reminiscent of unimodular gravity [66] . 
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